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ABSTRACT 

We discuss the polarization selection rules for splitting of the two principal 
electromagnetic modes that propagate in a vacuum polarized by a superstrong 
magnetic field (B > 0.1S cr ~ 4 x 10 12 G). We show that below the threshold 
of free pair creation the selection rules found by Adler in the limit of weak 
dispersion remain unaffected by taking the resonant effects into consideration, 
i.e., splitting of one mode is strictly forbidden, while splitting of the other is 
allowed. 
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1. Introduction 



Current models for the generation of coherent radio emission of pulsars require 
formation of free e + e~ pairs in the pulsar magnetospheres (e.g., Ruderman & Sutherland 



1975|; Michel |1991| ; Melrose |1993| ; Usov |1996|) . Such pairs may be created via conversion 



of 7-rays into e + e~ pairs in a strong magnetic field, ^ + B^e + + e~ + B (Erber |1966| ). 
Besides, some (if not all) pulsars are powerful sources of 7-ray emission (e.g., Hartmann 
1995 ; Thompson et al. |1995 ; Daugherty & Harding 1996 ). Therefore, propagation of 



7-rays in the magnetosphere is one of the central problems in the theory of pulsars (Usov 
& Melrose |I995| , [1996| ; Harding, Baring & Gonthier [IW7| ; Baring & Harding |200ll and 
references therein). 

In the magnetospheres of pulsars the plasma dispersion for 7-rays is negligible, and to 
understand the process of 7-ray propagation in the vicinity of a pulsar it suffices to consider 



propagation in the vacuum polarized by a strong magnetic field (e.g. Usov & Melrose |1995 
Bulik |1998| ). The principal modes of propagation for a photon in the magnetized vacuum 
are linearly polarized with the photon electric field either perpendicular (_L mode) or 
parallel (|| mode) to the plane formed by the external magnetic field B and the photon wave 
vector K. Both modes are generated in the magnetospheres of pulsars (Baring & Harding 



2001). A single photon in either mode can decay into a free e + e pair provided that a 



relevant threshold, e 1 sin^ = 2mc 2 for || mode and e 7 sin^ = mc 2 r] for _L mode, is exceeded, 
where e 7 is the photon energy, B cr = m 2 c 3 /eh = 4.4 x 10 13 G, r] = 1 + [1 + (2B/B CI )} 1 ^ 2 , 
and •§ is the angle between the photon wave vector K and the magnetic field B. For 
typical magnetic fields of pulsars, B ~ 10 12 G, the absorption coefficient for photo-pair 
production is at least a few orders of magnitude higher than the absorption coefficient 
for any other inelastic interaction of photons and magnetic fields (Adler |1971| ). However, 
for all known pulsars, 7-rays generated near the neutron star surface are produced in a 
state below the pair creation threshold (Usov & Melrose pL995| ; Baring & Harding |2001|) . 
In this case, the process of pair creation by single photons is kinematically forbidden, and 
the main (inelastic) interaction of photons and the magnetic field is a splitting into two 
photons 7 + 5^ 7l +72 + 5 (Adler |l97l| ; Stoneham [1979| ; Melrose |l9g3| ; Baring [19911) . 



Adler (|1971|) showed that in the limit of weak dispersion, where the refractive indices for 
both modes are very close to unity, only _L mode may undergo the splitting process. Near 
the threshold of pair creation, however, the dispersion law differs considerably from the 
vacuum case, |K| = u/c, and the weak- dispersion limit is no longer applicable (Shabad 
1975| ). The polarization selection rules for photon splitting with taking into account the 
resonant dispersion law near the threshold of free pair creation were discussed by Usov & 
Shabad ( |1983| ). Later, it was shown that bound e + e~ pairs (positronium) may be created 
by single photons in a strong magnetic field (Shabad & Usov |1985| , |1986| ; Herold, Ruder, & 
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Wunner [1985| ; Usov & Shabad [19851 ; Usov & Melrose fL99§ pJ96[) . The threshold of bound 
pair creation is slightly below the threshold of free pair creation, and therefore, creation of 
bound pairs may be allowed while creation of free pairs is forbidden yet. The dispersion 
curves corrected by taking into account formation of the mixed photon-positronium states 
were calculated in (Shabad & Usov |1986| ). These dispersion curves near and below the 
thresholds of pair creation differ significantly from the photon dispersion curves calculated 
by Shabad (|1975|) and used in (Usov & Shabad |1983j ). In this paper, using the corrected 
dispersion curves (Shabad & Usov |1986|) the polarization selection rules for splitting of both 
modes in a superstrong (B > 0.15 cr ) magnetic field are considered. 



2. The polarization selection rules for splitting 

In a time-independent, spatially uniform external magnetic field B, the dispersion 
curve for the principal mode of polarization a (a = {||, _L}) may be written in the following 
general form (e.g., Shabad |1975| ; Usov & Shabad |1983|) 



(u/c) 2 - K\ = f a {K\ 



where u = e 1 j% is the frequency of state (photon, positronium atom or their mixture), K\\ 
and K± are the components of the wave vector along and across the magnetic field B, and 
f a (K±) is a certain function of K\. In a superstrong magnetic field for the || mode below 
the threshold of free pair creation, {uj/c) 2 — K 2 < (2mc 2 /h) 2 , we have (Usov & Shabad 
1985] Shabad & Usov |T985l, Pgg) 
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where 



A(Kl) = ^oo(0,Kl)|^(0)| 2 exp (~^jp) , (3) 

£ O o(0, K±) = 2mc 2 — Ae o(0, Kj_) is the rest energy of positronium in the ground state when 
the Coulomb quantum number and the Landau numbers of the electron and the positron 
are zero (n c = = n + = 0), 
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Ae O o(0, K]_) = a 2 mc 2 



In 



r 



l + r 2 K 2 



(4) 



is the binding energy, V(0) is the longitudinal wave function of the positronium with the 
constituent particles taken at the same points (z + = z~), 

l^(0)| 2 = -m r a == , (5) 
a r^l + rlKl 

a = h/me 2 is the Bohr radius, r B = (ch/eB) 1 ^ 2 is the radius of the electron orbit, and 
a = e 2 /ch is the fine structure constant. The dispersion curve of the || mode which 
is given by equations (l)-(5) is photon-like for K\ < (2mc/h) 2 , and positronium-like 
for K\ > (2mc/h) 2 (see Fig. 1). At K\ ~ (2mc/h) 2 , this curve describes the mixed 
photon-positronium state. 

The positronium states of the lowest series (n + = n_ = and n c > 0), including the 
boundary of its continuum {uj/c) 2 — K 2 = (2mc/h) 2 , do not contribute into the polarization 
operator of a _L-polarized state (Shabad & Usov |1986"1) . Therefore, the dispersion curve 



(uj/c) 2 — Km = f±{Kf_) for the _L mode crosses the positronium spectra of the lowest series 
without interfering with them (see Fig. 1). As K\ grows further, this curve approaches 
eventually the dispersion curve of positronium at the state with n + + n_ = 1 and n c = 
when either the electron or the positron is at the first excited Landau level (Shabad & Usov 



19861) . 

At B > 0.15 cr , the geometrical optics or adiabatic approximation is valid for 



propagation of 7-rays in the pulsar magnetospheres (Usov & Melrose |1995| and references 
therein). In this case, the curvature of the magnetic field lines causes 7-rays of both modes 
to shift along the relevant dispersion curves (see Fig. 1). In the process of propagation, 
7-rays that are created significantly below the threshold of pair creation can pass through 
the mixed photon-positronium states and turns into positronium at \K±\ > 2mc/h. In 
this section, we discuss the polarization selection rules for splitting of both polarization 
modes irrespective of that the the initial state is a photon or a positronium atom. We 
consider the case when the initial state is below the threshold of free pair creation, 
(u/c) 2 - K 2 < (2mc 2 /h) 2 . 

The splitting reactions that may be in principle in a superstrong magnetic field are 

|| — > || + ||, |HI + - 1 ) IH-L + ll> IH-L + -L> -L-^-L + -U + ll> -L-HI + x and 

-L— >|| + || . At first we consider the kinematics of the splitting reaction ||— >\\ + || where the 
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initial and final states are the || mode. We use the frame where K» = for the initial state. 
If the energy of the initial state is hu < £oo(0,4m 2 c 2 ) this reaction is the photon splitting 
7 + B — > 7x + 72 + B, while at £oo(0, 4m 2 c 2 ) < Hu < 2mc 2 it is the positronium annihilation 
(e + , e~) + B — > 71 + 72 + B (see Fig. 1). Below, both these processes are called splitting 
into photons. 

The splitting reaction into two photons is kinematically allowed only if energy- 
momentum conservation can be satisfied: 

uj = uji + u>2 and K = K x + K 2 , (6) 

In our case, when K» = 0, we have K\\\ = —K 2 \\- 

Using equations (2) - (5), we can solve the dispersion relation (uj/c) 2 — K 2 = f\\(K]_) 



and find \K±\ as a function of u — ^Jlu 2 — c 2 K 2 , \K±\ = <p\\(u>), where u is the frequency 
in the frame where Kn is equal to zero. The function ip\\(&) is equal to zero at Q = and 



increases monotonously with increase of u (Shabad & Usov |1986|) . From equations (2) - (5), 
it follows that the deviation of \K±\ from the vacuum dispersion curve \K±\ = (uj/c) 
to the side of large values of \K±\ increases with increase of Co throughout the interval 
< Co < 2mc 2 /h (see Fig.l). In other words, we have <p\\(uji) > ^11(^2) for Q\ > 0J2 and 
<P\\(uii + uj 2 ) > <f\\(tii) + ^11(^2) for any values of and 0j 2 at Q\ + C02 < 2mc 2 /%. These 
properties of <p\\{£j) and equation (6) yield 



\K\ I = 



<p\\{u) = <p\\{ui + U2) > <f\\(uJi) + (p\\{u} 2 ) 



> n (Jul - c 2 K 2 ) + n (Ju 2 2 - c 2 K? 



= (p\\(u>i) + <p\\(&2) = \Ki,±\ + \K 2 ,±\ ■ (7) 

Therefore, the conservation of K± cannot be satisfied for the splitting reaction ||— >|| + || 
below the threshold of free pair creation, and this reaction is kinematically forbidden. This 
is valid irrespective of that the initial state is either a photon or a positronium atom. 

The splitting reactions ||— >|| + _L, ||— >J- + || and ||— >J_ + J_ are also kinematically 
forbidden below the threshold of free pair creation. Indeed, if we rewrite the dispersion 
relation {uj/c) 2 — K 2 = f±(Kj_) for the _L mode in the form \K±\ = (p±(u), the inequality 
<P\\(uj) > <£_!_(&) holds throughout the interval < to < 2mc 2 /h (Shabad [L972| , |1975| Shabad 
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& Usov 1986 ). Using this inequality the chain of inequalities (7) may be continued for 
splitting reactions when the initial state is the || mode and at least one final state is the 
_L mode. Finally, we have \K±\ > \Kx t ±\ + \K 2j ± | for all these reactions. Hence, splitting 
of the || mode is strictly forbidden below the threshold of free pair creation, i.e., Adler's 
conclusion remains unaffected by taking the resonant effects into consideration. 

Since the dispertion curve \K±\ = <p±(u)) for the _L mode crosses the positronium 
spectra of the lowest series (n + = n_ = and n c > 0) without interfering with them, the 
resonant effects are unessential for splitting of the _L mode below the threshold of free 
pair creation, Co < 2mc 2 /Ti (see Fig. 1). In this case, the splitting reaction _L— »_L + _L is 
kinematically forbidden while the splitting reactions J_— >|| + ||, J_— >|| + _L and _L— >_L + || 
are kinematically allowed (Adler |1971| ). Anyone can also come to this conclusion from 
the following properties of the functions (p±(&) and <p\\{&) below the threshold of free 
pair creation, cp±(Qi) > <px.{u>j) for Q\ > cD 2 , <fx(u>i + 0J2) > <P±(&i) + <P±(&2), and 

<P\\(u>) > (p±(&) (Shabad & Usov |1986|) . The splitting reactions I — ^ | ] + J_ and J_— >± + || 

that involve an odd number of ||-polarized photons are forbidden by CP invariance in 
the limit of zero dispersion (Adler 1971). In the magnetized vacuum, dispersive effects 



guarantee a nonzero probability for these reactions. Therefore, the channels _L— >\\ + _L 
and _L— >_L + || are strongly suppressed in comparison with T— >|| + || when the refraction 
index of the magnetized vacuum is close to unity, n — 1 ~ 0.1a(B/B CI ) 2 1 (Adler 
|1971[ ). However, this suppression may not hold in an extremely strong magnetic fields 
(B ^ (10/a) 1 / 2 ^ ~ 10 15 G). 



3. Discussion 

In this Letter, we have discussed splitting of the principal modes that propagate in 
the vacuum polarized by a superstrong (B > 0.1S cr ) magnetic field. We have shown that 
splitting of the || mode is strictly forbidden below the threshold of free pair creation. 
This is valid irrespective of that the initial state is either a photon (\K±\ < 2mc/h) or a 
positronium atom (\K±\ > ImcjK). Splitting of the _L mode is allowed, and only the one 
channel (-L— >|| + ||) of splitting operates if the magnetic field is not extremely high (Adler 



19711) . 



In the study of the polarization selection rules for splitting, Adler ( |1971| ) has used the 
refractive indices 
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_ C |K|_ yj, ± (o>) + Kf 

while we use the dispersion curves = tp\\ t ±_{uj). From equation (8), the properties of 
nn x used by Adler ( |1971| ) follow immediately from the properties of (pu± used in our study, 
and therefore, both these methods are equivalent. 

The fact that the resonant effects do not change the polarization selection rules found 
by Adler for photon splitting in the limit of weak dispersion may be easily explained. 
Indeed, to formulate the selection rules for the case when photons are significantly below 
the threshold of pair creation the following properties of functions (p\\(u)) and <p±(ui) has 
been used, in fact, by Adler ( |1971| ), <p\\(u>) > <p±(u) and for both principal modes the 
deviation of \K±\ from the vacuum dispersion curve = (£/c) to the side of large values 
of \K±\ increases with increase of to. These properties only strengthen by the resonant 
effects, especially when the dispertion curve of the || mode approaches the positronium 
state at \K±\ ^> 2mc/h (see Fig. 1). Therefore, for a positronium atom at the ground 
state (n + = n_ = n c = 0) with ^> 2mc/h its stability against spontaneous decay into 
photons is evident, and it has been mentioned long ago (Herold et al. |1985| ). 



The dispersion curves calculated in (Shabad & Usov |1986|) and used in this Letter 
are exact relativistic consequences of the adiabatic approximation, and therefore, they are 
valid for an arbitrary strong magnetic field. [This relates to equations (1), (2) and (3).] 
The nonrelativity is introduced only when the binding energy of the positronium (4) and 
the longitudinal wave function (5) that are the solutions of the nonrelativistic Schrodinger 
equation are substituted into equation (3). The relativistic corrections to both the binding 
energy and the longitudinal wave function are small if the positron and the electron remain 
nonrelativistic along the field direction, i.e., if Ae o(0,0) <C mc 2 (Angelie & Deutsch |1978| ; 



Lai & Salpeter 1995 ). From equation (4), it follows that the dispersion curves calculated in 



(Shabad & Usov |1986|) and therefore, our consideration are valid even if the magnetic field 
is as high as the virial magnetic field (~ 10 17 — 10 18 G) for neutron stars. 

There is now compelling evidence for the existence of "magnetars" - neutron stars 
with magnetic field strenghths in excess of B cr . The evidence comes primarily from 
observations of rapid spin down in the pulsations observed from soft 7-ray repeaters (e.g., 
Kouveliotou et al. |1998j , |1999| ). Recently, radio surveys have discovered a few pulsars 



with the surface magnetic fields approaching 10 14 (Camilo et al. 2000). Kulkarni & Frail 
( |1993|) and Kouveliotou et al. ( |1994j ) estimate that the birthrate of magnetars is about 
10% that of ordinary pulsars. The processes of photon splitting and bound pair creation 
may be very important for magnetars (Thompson & Duncan |1995| ; Usov & Melrose |1995 
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Baring & Harding 2001 and references therein). We briefly discuss these processes in the 



magnetospheres of pulsars with very strong magnetic fields. 

For pulsars with the surface magnetic field B$ ^ 10 13 G, _L-polarized 7-rays generated 
near the neutron star surface split into || -polarized 7-rays before they reach the pair creation 
threshold (e.g., Baring & Harding [200 1|) . In turn, in such a strong magnetic field ||-polarized 
7-rays create bound pairs which are stable, in the absence of such external factors as electric 



fields and ionizing radiation (Usov & Melrose |1995| and references therein). The bound pairs 



form a gas of electroneutral particles. Such a gas does not undergo plasma processes, like 
plasma instabilities, which are responsible for the generation of coherent radio emission of 
pulsars. Maybe, the suppression of free pair creation in superstrong magnetic fields results 
in a death line of pulsars at B$ ~ B cr (e.g., Baring & Harding |1998| ; Heyl & Kulkarni |1998| ). 



For this suppression it is very important that splitting of _L-polarized 7-rays into || -polarized 
7-rays prevents from formation of bound pairs at the state with n + + n_ = 1 when either 
the electron or the positron is in an excited state with the spin quantum number, s — 1, 
opposite to that in the ground state. Otherwise, bound pairs in the excited state may be 
ionized by the recoil from photons radiated in a spin-flip transition, and free pairs may form 
in the magnetospheres of pulsars. 

Since bound pairs do not screen the electric field E\\ , the formation of bound rather 
than free pairs may cause increases in the height of the polar gap and the total power 
carried away by relativistic particles from the polar gap into the pulsar magnetosphere by 
large factors (e.g., Shabad & Usov p~985|) . For the case when there is no emission of particles 



from the surface of a strongly magnetized neutron star, these increases were considered by 
Usov & Melrose ( |1995| , |1996| ) in detail. Recently, the polar gap model, which is rather close 
to reality for typical pulsars, is developed for the case when particles flow freely from the 
stellar surface (Muslimov & Harding |1997| ; Harding Sz Muslimov |1998| ; Zhang & Harding 
2000| ; Hibschman & Arons |2001| ). In this model, the E\\ field in the polar gap is due to the 



effect of general relativistic frame dragging discovered by Muslimov & Tsygan (|1992|). For 



the pulsar period P ^ 0.1 s and 0.15 cr <^ Bo, 10-B cr , this field (e.g., Harding & Muslimov 
1998|) is significantly smaller that the value E\? a ~(1 — 2)xl0 10 V cm -1 , at which filed 



ionization of bound pairs becomes important (Usov & Melrose |1996|) . The mean free path 
of bound pair photoionization is Z p h ~ 10 5 (r/10 2 ) 3 (Tg/10 6 K) -2 cm (Usov & Melrose 



1995| ), where V is the Lorentz factor of bound pairs and Tg is the surface temperature. 
For high-energy bound pairs (r > 10 3 — 10 4 ) that might be responsible for the E'li-field 
screening after their ionization, we have that is larger than the polar gap height even 
if Ts is as high as a few x 10 6 K. Therefore, the polar gap model developed in (Muslimov 



& Harding |1997| ; Harding & Muslimov |1998| ; Zhang & Harding |2000] ; Hibschman & Arons 



2001|) is probably accurate for magnetic fields B$ ^ (0.1 — 0.2)5 cr while for B$ > 0.2B cr it 
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is necessary to modify this model by taking into account photon splitting and bound pair 
creation. 

In the vicinity of pulsars with B$ ^ 10 13 G splitting of _L-polarized 7-rays into 
1 1 -polarized 7-rays provides a mechanism for the production of linearly polarized 7-rays. 
At energies e 7 <; 10 2 MeV, the 7-ray polarization may be up to 100% (Shabad & Usov 



1983| ). By observing the polarization of the 7-ray emission of pulsars it would be possible 



to estimate the magnetic field near the pulsar surface. 
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Fig. 1. — Dispersion curves for the || mode (thick solid line) and for the _L mode (thick dashed 
line). The sloped straight line is the vacuum dispersion curve, |Jfj_| = uj/c. Positronium 
dispersion curves for the states with n+ = n_ = n c = and with n + + n_ = 1 and n c = 
are shown by the dotted line and by dot-dashed line, respectively. The horizontal straight 
dashed lines are the boundaries of continua for the first two series of bound states with 
n_|_ = n_ = or n + + n_ = 1. There is an infinite number of branches for the states with 
n c > 1 between the continuum boundaries and the positronium dispertion curves for the 
states with n c = (Shabad & Usov |1985| , |1986|) . They are not shown for simplicity. 



